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Abstract 

C/^ \ New spectroscopic observations of the double-lined eclipsing binary AQ Cas 

^' 

O ! been analysed for the fundamental properties of the components. Analyses 

ii3 ■ 

t^ ' show that the system consists of a massive primary with a mass of 17.63±0.91 

Mq and radius of 13.48±O.64R0 and a secondary with 12.56±0.81 Mq and 

>• . radius of 23.55±0.73 R0, corresponding spectral types of B0.5(±2) II-III + 

(N 

0\ \ B3(±l) II. The secondary star fills its corresponding Roche lobe and mass 

'nI" , 

^ ■ transfer to the primary star is going on. This stream considerably does affect 

^' 

O ■ the photometric observations both starting from the second quarter up to 

cn 

the first contact of primary eclipse and just at the second maximum. Thus, 
the light curve is distorted and tightly depended on the wavelength of the 



H ■ observations. The available multi passband light curves have been analysed 

by taking the stream effects, as either hot or cool spots, into account. The 
comparison of the models and observations in the log (L/Lq) - log Tg// 
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and log g — log Tf.ff diagrams clearly shows that the more massive star is 

consistent with models and is predicted to be close to the phase of hydrogen 

shell ignition. Average distance to the system is estimated as 4150±240 pc 

using the BVJHK magnitudes and V-passband extinction. 

Keywords: Binaries Eclipsing - stars: fundamental parameters Individual 

method:spectroscopy 

1. Introduction 

A massive star with a mass higher than 8-9 M0 evolves through 
all cycles of nuclear burning up to iron when a core collapse su- 
pernova is formed. Such high mass stars undergo a strong loss of 
mass by the stellar wind which decreases their mass during nuclear evo- 
lution. Observations of many massive stars have shown that the mass loss 
is larger for the more massive stars. The mass losses are so large that the 
outer atmospheres are removed and the convective cores may become visible 
at the surface. In this way some elements processed by nuclear reactions are 
carried out to the surface of a star. Spectroscopic observations of massive 
stars at advanced evolutionary stages reveal that the outer layers are en- 
hanced in helium, carbon, nitrogen and oxygen. The evolutionary tracks are 
considerably different for the mass losing stars than those at constant mass 



( jDe Loore fc Dooml . Il992l ). When a star has a close companion its structure 
and evolution are significantly changed. While mass and chemical compo- 
sition of a star define its evolution in the case of a binary system the total 
mass, the mass-ratio and the orbital period would affect its evolution. The 
mutual interaction between the components may cause to mass and angu- 



lar momentum exchange which result in changes in the mass-ratio and the 
orbital period of the system. Although a number of the differences between 
the single and binary stars are well understood several are still in discussion, 
esp ecially for massive binary s t ars. 



Podsiadlowski. Joss &: Hsul ( Il992l ) have investigated in detail how binary 
interaction affects the pre-supernova evolution of massive close binaries and 
the resulting supernova explosions. About 30 percent of all massive stars 
are substantially affected by interactions with close binary companions. 
Further evolution of a close bi nary system is ti ghtly depended upon the mass 



and angular momentum loss. 



EggletonI ( I2OOOI ) called attention, for the first 



time, that mass ar id sub s equeii t angular momentum loss significantly affect 



binary evolution. 



Sarna] ( 119931 ) have shown that during evolution of /3 Per 



15% of the initial total mass about two of third of the mass lost and 30 % 
of the total angu l ar mo mentum were lost. Van Rensbergen his collaborators 
van Rensbergeru (120031 ) have shown that conservative calculations cannot 
produce Alg ols with large mass-ratios. 



Recently 



Torres. Andersen. &: GimenezI (120 lOf ) collected accurate masses 



and radii for the components of 95 double-lined detached binary systems. 
Their criterion was such that the mass and radius of both stars should be 
determined within an accuracy of 3 percent or better. The number of stars 
is decreasing in both directions towards the massive stars (M>10 M©) and 
low-mass stars (M<1 Mq). In the last three years the number of low-mass 
stars for which mass and radius deter r ained a re considerably increased(see for 
example ( ICakirh. Ibanoglu.fc Sipahil . l2013al Jbl). Still there are some debates 
about the formation mechanism and evolution of massive binaries. Therefore, 



we initiated a spectroscopic study for the massive binary systems. 

AQ Cas is a semi-detached, /3 Lyrae type echpsing binary with an orbital 



Struvd ( Il946l ) and the 



period of 11.721 days. Its spectra were obtained by 

compone nts were cla ssifie d as B3 and B9. F irst photometric observations 

made by 



Olsonl f ll985h and 



Nha fc ChunI (119881 ). 



iater on UBV observations 



LeeetaL 



fll993h and 



Olson 



and I (Kron)uvby light curves were published by 
(119941 ). respectively. There is large light variations between eclipses and 
light curve has a deep totality lasting about ten hours. In addition, light 
curves are significantly asymmetric which appears to tightly depended on the 
wa yelength. About 40 years later some spectra of the system were obtained 



by 



LeeetaL 



(Il993[ ). These spectra showed the traces of both stream and 
rotation effects. 

We report here spectroscopic observations which are added to the pre- 
vious data and analysed together with the photometric observations for ob- 
taining accurate masses and radii of the components. In addition, we have 
interpreted all available photometric and spectroscopic data on AQ Cas for 
understanding its present structure and evolution as well as physical proper- 
ties of the stars. We also discuss the plausible causes of asymmetries in the 
observed multi passband light curves of the system. 



2. Spectroscopic observations 

Optical spectroscopic observations of the AQ Cas were obtained with the 
Turkish Faint Object Spectrograph Camera (TFOSCJj attached to the 1.5 



^http://tug. tug. tubitak.gov.tr/rttl50_tfosc.php 



m telescope in July, 2010 under good seeing conditions. Further details on the 



telescope and the spectrograph can be found at, http://www.tug.tubitak.gov.tr 



The wavelength coverage of each spectrum was 4000-9000 A in 12 orders, with 
a resolving power of A/AA ~7000 at 6563 A and an average signal-to-noise 
ratio (S/N) was ~120. We also obtained high S/N spectra of several early 
type standard stars 21Cyg, rHer, HR153 and 21 Peg for use as templates 
in derivation of the radial and rotational velocities. The electronic bias was 
removed from each image and we used the 'crreject' (IRAF-immatch- 
imcombine-"Type of rejection") combineoption for cosmic ray removal. 
Thus, the resulting spectra were largely cleaned from the cosmic rays. The 
echelle spectra were extracted and wavele ngth calibrated by usin g Fe-Ar lamp 



source with help of the IRAF echelle (JTonry k, David . Il979l ) package 



2.1. Radial velocities 

To derive the radial velocities, the eleven spectra obtained for the system 
are cross-correlated against the template spectra of st andard stars on an 



order- by-order basis using the FXCOR package in IRAF ( ISimkiru . Il974f ). 

The spectra showed two distinct cross-correlation peaks in the quadra- 
tures, one for each component of the binary. Thus, both peaks are fitted 
independently with a Gaussian profile to measure the velocities and their 
errors for the individual components. If the two peaks appear blended, a 
double Gaussian was applied to the combined profile using de-blend func- 
tion in the task. For each of the eleven observations we then determined a 
weighted-average radial velocity for each star from all orders without signifi- 
cant contamination by telluric absorption features. Here we used as weights 
the inverse of the variance of the radial velocity measurements in each order. 



as reported by FXCOR. 

Since the orbital period of the system is considerably longer and the 
totality lasts about ten hours it is hard to obtain reliable times of mid-eclipse. 
Therefore, only nine photoelectric times for primary eclipse are collected from 
the GATEWAY data-base. A linear least-squares fit to the data yields the 
following ephemeris 

Mini = HJD 2453920.897(27) + ll'^.72107(18) x E 

where the bracketed quantity is the uncertainty in the last digits of the 
preceding number. 

The heliocentric radial velocities for the primary (Vp) and the secondary 
[Vs) components are listed in Table 1 , along with the dates of observations 
and the corresponding orbital phases computed with the new ephemeris given 
above. The velocities in this table have been corrected to the heliocentric 
reference system by adopting a radial velocity value for the template stars. 
The radial velocities are plotted again st the orbi t al ph ase in Fig. 1, together 



with the radial velocities obtained by 



Leeetal 



(Il993[ ) (empty circles). The 



three velocities of the primary star in phases 0.12-0.16 and nine velocities 
in phases 0.81-0.90 are affected from gas stream and rotation. In addition, 
the velocities of both components in the orbital phase about 0.38 are sys- 
tematically deviated from the orbit. Therefore, thirteen measurements for 
the primary and one for the secondary are not plotted in Fig. 1. There is no 
systematic difference between our and their remaining measurements. 

We analysed all the radial velociti es for the initial orbital para meters 



using the RVSIM software program (JKane. Schneider, fc Gel . 120071 ). Fig- 



ure 1 shows the best-fit orbital solution to the radial velocity data. The re- 
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Figure 1: Radial velocities for the components of AQ Cas. Symbols with error bars show 



the radial velocity measurements for the components of t he system 



Lee et al. 



primary: filled circles, 
(|l993l ) are shown by 



secondary: open squares). The velocities measured by 
empty circles. 



suits of the analysis are as follow: 7= -16.6±2.2 km s~^, i^i=120.9±3.9 and 
i^2=169.6±3.4 km s~^with circular orbit. Using these values we estimate the 
projected orbital semi-major axis and mass ratio as: asini=67.4±1.2 Rq and 



|^=0.713±0.027. 

Ml 



2.2. Rotational velocity 

The width of the cross- correlation pr o file is a good tool for the mea- 
surement of vsini (see, e.g.. iQueloz et al.l (Il998l )). The rotational velocities 
(vsini) of the two components were obtained by measuring the FWHM of 
the CCF peaks in nine high-S/N spectra of AQ Cas acquired close to the 
quadratures, where the spectral lines have the largest Doppler-shifts. In or- 
der to construct a calibration curve FWHM-f sini, we have used average 
spectra of the l Psc and r Her. The FWHM of the CCF peak was measured 
and the FWHM-w sin i calibration was established. The v sin i values of the 



Table 1: Heliocentric radial velocities of AQ Cas. The 
columns give the heliocentric Julian date, the orbital 
phase (according to the ephemeris in §2), the radial ve- 
locities of the two components with the corresponding 
standard deviations. 

HJD 2400000+ Phase Star 1 Star 2 

Vri a Vs cr 



56131.4953 


0.6010 


77.5 


3.2 


-118.6 


4.3 


56131.5470 


0.6054 


73.2 


3.4 


-120.6 


4.6 


56133.5279 


0.7744 


121.4 


2.8 


-174.0 


3.8 


56134.5740 


0.8637 


90.1 


3.1 


-140.5 


4.2 


56137.5130 


0.1144 


-110.4 


3.1 


101.9 


4.4 


56326.5500 


0.2423 


-144.7 


2.1 


155.3 


3.0 



components of AQ Cas were derived from the FWHM of their CCF peaks. 
We find 287±5km s^^ for the primary and 98±9km s^^ for the secondary 
star. 

3. Light curves and their analyses 

Th e first uvhyl photoelectric observations of AQ Cas were obtained by 



Olson! (Il985[ ). The resulting light curves are significantly variable in the 
sense that the light just before the first contact is depressed by about 0.1 
mag in the u-passband relative to the fourth contact. This difference is de- 
creased towards longer wavelengths. The light loss before the first contact is 
attributed to the stream seen projected against the gainer. In addition, some 
light fluctuations are seen beginning fror a mid-secondary-ecli pse. Additional 



observations of AQ Cas were obtained by 



Olson &: Belli (119891 ) using the same 



filter set. First wide passband UBV observations obtained between 1982 and 



1989 were published by iLee et al.l (Il993l ). Despite the large fluctuations in 
fluxes the average shape of the light curve is revealed. Their BV light curves 
are almost symmetric except the secondary eclipse. The ingress is steeper 
than the egress. Five- color interme diate-band uvhyl observations of the sys- 
tem were obtained by lOlsonI (Il994| ). Thus, almost complete multi- passband 



light curves were revealed. The asymmetry in the light curves was increasing 
towards the shorter wavelengths. In addition, the asymmetry is in an oppo- 
site direction as we go to the longer wavelengths. In the u-passband light 
curve was the most asymmetrical shape, being the first maximum dimmer 
than the second. However, the I-passband light curve is almost symmetric, 
the light level at the first maximum is slightly higher than that of the second 
maximum. 

The observations of AQ Cas were also made by two familiar surveys. The 
R-passband l ight curve was obtain ed by the Northern Sky Variability Sur- 
vey (NSVS) (jWozniak et al.l . |2004J ). The light curve is very similar to that 
obtained by Olson in the /-passband. The last V- passband light curve of 
AQ Cas was obtained by the Int ernational Gamma- Ray Astr ophysics Lab- 



oratory {INTEGRAL) mission flAlfonso-Garzon et al.l . I2OI2I ). With 3139 
V-passband measurements the shape of the light variation is clearly defined. 
The magnitude of the system appears to reach maximum just at phase 0.75 
and a steady decline up to the first contact is seen. The egress of the primary 
eclipse is steeper than that of egress. In contrary the egress of the secondary 
eclipse is steeper than that of egress. Both the minima and maxima are 
asymmetric. In Fig.2 we show Olson-I Lee-V, NSVS-R and INTEGRAL-V 
light curves of AQ Cas. 




0.7 G.B 0.9 



Phase 



Figure 2: The Olson's /, Lee's V, NSVS-i? and INTEGRAL-V light curves of AQCas. 
The continuous lines show the best-fit model. 
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Since the analysis of the hght curves is depended on the effective tem- 
perature of the hotter component try to constrain the effective temperature 
and spectral type of the primary star using the uvhyl and JHK magnitudes. 
The interme diate passband colors of both components are already given by 
Olson! ( 119851 ) with an interstellar reddening of E(b-y)=0.64mag. V=10.31, 
B-V=0.53, J- H=0.182±0.039 and H-K=0.170±0.031 mag were taken from 
SIMBAD data-base and 2MASS catalog. Using the colors and reddening 
we estimate a spectral type and effective temperature of the primary star as 



B0.5II1 and Teg = 2 7 000±1000K bv color-spectr a 



relatio ns g iven by 



teoOOh and 



de 



Tokunaga 



ager 



mm 



fc Nieuwenhuijzen 



type -e ffective temperature 



fll987h 



Drilling &: Landolt 



We used the eclipsing binary light curve mod eling algorithm o f 



( Il97ll ). as implemented in the PHOEBE code of IPrsa k Zwitterl (120051 ). The 



Wilson fc Devinney 



non-linear sq uare- root limb-darkening and bolometric limb-darkeii i ng ar e 
adopted from iDiaz-Cordoves. Claret. fc GimenezI (119951 ) Ivan Hammd (119931 ) , 
respectively. The gravity-brightening coefficients gi=g2=l-^ and albedos 
Ai=A2=l.^ were assumed for the components. The rotational velocity of the 
primary star is measured about 287 km s^^, being five times faster than syn- 
chronous rotation. While the less massive star is appeared to be synchronous 
with the orbital one, for the more massive component -Fi=5 is assumed. Us- 
ing a trial-error method, we obtained a set of parameters which represented 
the observed light curves. This preliminary analysis shows that the sec- 
ondary component fills its corresponding Roche lobe. Therefore, Mode-5 
was adopted. 

The shapes of the light curves are tightly depended on the wavelength. 
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Photometric and spectroscopic observations clearly show that there is a mass 
stream from the secondary to the primary star. The stream impacts directly 
to the more massive star. This impact region is seen by the observer just 
at the second quarter. Therefore, additional light is observed. Thereafter, 
the stream begins to absorb the light of the primary star which causes light 
loss. On the other hand, during the second quarter the observer receives 
light from the steam in addition to both stars. Therefore, the brightness of 
the system is slightly increased. All these effects can only be modelled by 
the hot or cool spots in the Wilson-Devinney code. 

The adjustable parameters in the light curves fitting were the orbital in- 
clination, the surface potential of the primary star, the effective temperature 
of the secondary, and the wavelength-dependent luminosity of the hotter star, 
the zero-epoch offset. Our final results are listed in Table 2. The uncertainties 
assigned to the adjusted parameters are the internal errors provided directly 
by the code. The computed light curve corresponding to the simultaneous 
light- velocity solution is compared with the observations in Fig. 2. 

Combining the results of radial velocities and light curves analyses we 
have calculated the absolute parameters of the stars. The fundamental stel- 
lar parameters for the components such as masses, radii, luminosities are 
listed in Table 3 together with their formal standard deviations. The stan- 
dard deviations of the parameters have been determined by JKTABSDIMJf 
code, which calculates distance and other p hysical parameters using several 



different sources of bolometric corrections (ISouthworth et al. 



20051). The 



This can be obtained from http://http://www. astro. keele.ac.uk/~jkt/codes. html 
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Table 2: Results of the simultaneous analyses of the Olson's /, Lee's V, NSVS-i? and 
INTEGRAL-V light curves for AQ Gas light curves. 



Parameters 



T./A (K) 
Te/A (K) 

r2 



(Olson-/) 



Adopted 



(^i^^%^f {INTEGRAL-V) 0.4424±0.0023 
a 0.032 



84.41±0.40 

17 000 [Fix] 

16 700±400 

5.781±0.166 

4.640±0.155 

0.1995±0.0088 

0.3484±0.0090 

0.4028±0.0101 

0.4584±0.0114 

0.4190±0.0119 



mass for the primary of Mp = 17.63 ± 0.91M q and secondary of Ms = 
12.56 ± O. SlMr;^ are consisting of a B 0.5II-III star and an evolved B3 super- 
giant star ( iDrilling fc Landolt 1 . 12000[ ) . 

The interstellar reddening is estimated by lOlsonI (jl985[ ) as E(b-y)=0.64 
mag which corresponds to E(B-V)=0.88 mag. Using th is value and the BVR- 



JHK magnitudes and bolometric corrections given by 



Girardi et al. 



(12002[ ) 



we estimated an average distance to the system as 4150±240 pc. 



4. Results and Conclusions 



First spectroscopic observations were made by IStruvd ( 119461 ) who esti- 



mated minimum mass for the components as llMfr^- f- 



Lee et al. 



5M(Ti and spectral type 



(Il993[ ) derived absolute 



of B3+B9. Almost a half century later 
physical properties of the component stars using their spectroscopic and pho- 
tometric observations. The masses and radii are found as I8.IM0+I2.9M0 
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Table 3: Properties of the AQ Cas components 



Parameter 



Primary 



Secondary 



Mass (Mq) 
Radius (Rq) 

Teff (K) 

log (L/Lq) 

log 9 {cgs) 

Spectral Type 

{vsin i)obs. (km s"-"^) 

{vsin i)calc. (km s'^) 

a(R0) 

V-y (km s-l) 

1 

d (pc) 



17.63±0.91 12.56±0.81 

13.48±0.64 23.55±0.73 

27 000±1000 16 700±400 

4.940±0.076 4.590±0.050 

3.425±0.039 2.793±0.026 

B0.5(±2)II-III B3(±1)II 

287±5 98±9 

15±1 4±1 

67.59±1.21 

-9±1 

0.713±0.027 

4150±240 



and I3.5R0+25.8R0, respectively. They classify the system as a semi-detached 
eclipsing binary with out any concl usion about the spectral types of the com- 



ponents. A year later 



m 



Olsonl (11994J ) published the five- color light curves and 



their analysis. He also measured radial velocities of the cooler star from the 
spectroscopic observations of 01 A7774 line. He estimates m asses for the com- 



ponents I6.6M0 and 4.72M0 , radii of 7.55R0 and I6.6R0. IPolushinal (120041 1 



published a catalog of massive close binaries with early-type components 
which includes AQ Cas. The spectral types and masses of the components are 
given as 08.5III+08.5IH and 29.5M0+24.6M0. Our re-analyses of all avail- 
able photometric and spectroscopic data yield the mass, absolute radius and 
effective temperature for the primary and secondary stars are 17.63±O.91M0, 
13.48±O.64R0 , 27000±1000, and 12.56±O.81M0, 23.55±0.73, 16 700±400 , 
respectively. These values correspond to the spectral type of B0.5IH and 
B3II. Thus, the masses of the components have been derived with an accu- 
racy of 5 percent for the primary and 6 percent for the secondary and radii 
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Figure 3: log {L/LQ)-\og T^ff (left panel) and log q — log Tf-jf (rig ht panel) plots of the 



primary component of AQ Cas compared with 



Ekstrom et al 



(12012) models for Z=0.014. 



of 5 percent and 3 percent, respectively. 

In Fig. 3 we plot the location of the primary star of AQ Cas in log (L/Lq)- 
log Teff (left panel) and log g — log T^jf (right panel) diagrams with er- 
ror bars. The e volut ionary tracks of a 17.63 Mq have been constructed by 



Ekstrom et al. 



mm models for X=0.720, Y=0.266 and Z=0.014. Stellar 
rotation is taken into account in their models. Despite their models provide 
a good description of the average evolution of non-interacting single stars, 
locations of the gainer are in good agreement with evolutionary tracks as 
seen in both panels. The primary star is predicted to be close to the phase 
of hydrogen shell ignition. 

AQ Cas is a classical Algol. While the more massive star is filling its Roche 
lobe by about 0.52 the less massive star completely fills and overflows its lobe. 
A mass transfer is taking place from the less massive (donor) star to the more 
massive (gainer). Mass flows from the inner Lagrangian point and the trajec- 
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tory and prop e rties of the mass-transferring stream are studied in detail by 
Lubow &: Shul (119751 ) . The type of accretion structures can be predicted from 
position of a gainer in the well-known radius — massratio diagram. If the 
orbital period of a system is too long, the fractional radius will be too small 
that mass transfer results in classical accretion disks. If the fractional ra- 
dius of the gainer is sufficiently large, stream directly strikes its photosphere. 
This impact leads to the formation of a hot impact region and variable ac- 
cretion structures. The gainer of the AQ Cas, with r=0.2 and q=0.7 , locat es 
in the region of impactors (see Fig. 3 in lPervisoglu. Tout.fc Ibanoglul ( I2OIOI )). 
Although its position is close to U Cep in the r — q plane, the mass-ratio 
and orbital period for AQ Cas are very different. It seems that AQ Cas is a 
mass-transfe rring system with higher mass - ratio and longer orbital period. 

Recently iDervisoglu. Tout.fc Ibanoglul (I2OIOI ) studied spin angular mo- 
mentum evolution of the accreting components of Algol-type interacting bi- 
naries. They demonstrate that a small amount of mass transfer leads to 
accumulation of enough angular momentum which spins the gainer up to 
its critical rotation velocity. Since the gainer is in Algols have radiative en- 
velopes, angular momentum and mass loss are not taken place due to the 
enhanced magnetic activity, as are occured in solar-type chromospherically 
active stars. Therefore, they take into account generation of magnetic fields 
in the radiative atmospheres in differentially rotating star and the possibil- 
ity of angular momentum loss driven by strong stellar winds. Differential 
rotation induced by the accretion itself may produce stellar winds which 
carry away enough angular momentum to reduce their rotational velocities 
to the presently observed values. They suggest that this self consistent model 
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should be more effective in the angular momentum loss of long-period Al- 
gols. The rota t ional velocity of the gainer was measured as 240 km s^^ by 



Etzel fc Olson 



m 



19931 ). and they estimate an Fi value of 7.3. We measured a 



rotational velocity of Vrotsin z=287±5km s"^ which is five times faster than 
synchronous rotation. As the mass of the gainer increases the rotation veloc- 
ity is decreasing, as shown in their Fig.2. The Fi value found in this study is 
in agreement with their estimate fo r massive primaries w ith orbital periods 



longer than 5 days. As discussed by llbanoglu et al.l (120061 ) interacting Algols 



with g > 0.3 and P > 5 days have almost the same angular momentum to 
those of detached binaries. So, with a mass-ratio of 0.7 and an orbital pe- 
riod of 11.7 days the angular momentum loss during the evolution of AQ Gas 
is not too high. Comparison with the evolutionary tracks indicates similar 
results. 
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